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THE  RECURRENCE  TENDENCY  AND  FORECASTS  OF  MAGNETIC  ACTIVITY 

A.  H.  Shapley 

Abstract- -The  recurrence  tendency  in  magnetic  character  figures  is  evaluated  by 
auto-correlation  coefficients  for  various  intervals  up  to  83  days  for  each  year  from  1890 
to  1944.  The  strength  of  the  27-day  recurrence  pulse  varies  during  the  55  years,  being 
generally  strong  from  about  two  years  after  sunspot  maximum  throi^h  the  minimum 
year.  Evidence  for  persistent  recurrence  intervals  other  than  multiples  of  about  27  days 
is  very  weak.  The  magnitudes  of  the  correlation  coefficients  for  short  intervals  indicate 
that  the  time  unit  of  geomagnetic  activity  is  two  to  three  days.  Formulas  for  the  predic¬ 
tion  of  magnetic  character  figures  by  the  method  of  least  squares  are  based  on  the  mean 
auto -correlation  coefficients  for  intervals  of  1>  26,  27,  28,  and  54  days  for  stages  of  the 
sunspot  cycle  with  weak  and  strong  recurrence  conditions.  A  prediction  formula  with  a 
non-linear  term  is  also  derived. 

The  recurrence  tendency- -The  tendency  for  magnetically  disturbed  or  quiet  periods  to  repeat 
themselves  at  intervals  of  about  27,  54  or  higher  multiples  of  27  days  is  not  equally  pronounced 
during  the  various  stages  of  the  solar  cycle.  Although  CHREE  and  STAGG  [see  *  References”  at 
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end  of  paper,  1927]  found  the  tendency  present  in  each  of  tiie  years  1906  to  1925,  the  recurrence 
was  better  defined  in  the  years  of  low  simspot  activity  and  also  in  the  years  when  the  average  sun- 
spot  latitude  was  low  [see  also  ARCHENHOLD,  1939].  BARTELS  [1932]  points  out  the  variabUity 
in  the  strength  of  the  tendency  as  shown  in  his  charts  devised  for  this  purpose.  MAUNDER  [1905], 
CHREE  and  STAGG  [1927],  and  GREAVES  and  NEWTON  [1929]  all  found  that  the  recurrence  ten¬ 
dency  was  stronger  for  the  moderate  types  of  magnetic  disturbances  than  for  the  great  storms. 
These  restdts  come  from  analyses  by  the  superposed  epoch  method  of  CHREE  [1912],  in  which  by 
selecting  disturbed  (or  quiet)  days  according  to  some  criterion  and  computing  the  average  mag¬ 
netic  activity  +  n  days  from  the  selected  day  it  is  possible  to  show  a  recurrence  tendency  without 
involving  a  preconception  of  the  length  of  the  interval.  Maunder,  and  Greaves  and  Newton  used  the 
amplitude  of  disturbance  as  the  basis  for  selection  of  days,  thereby  including  many  data  in  years 
when  many  disturbances  occurred  but  relatively  few  data  in  other  years.  Chree  and  Stagg,  and 
Archenhold  selected  the  five  most  disturbed  days  in  each  month,  thus  equalizing  the  data  from 
each  year,  but  varying  the  degree  of  storminess  of  the  days  selected.  The  former  also  made 
separate  use  of  the  five  quietest  days  per  month,  the  combination  giving  perhaps  a  more  uniform 
body  of  data.  The  selection  of  days  of  comparable  disturbance  on  the  one  hand  gives  an  unequal 
distribution  of  selected  days  during  the  various  phases  of  the  cycle  of  solar  or  magnetic  activity; 
on  the  other  hand  the  selection  of  a  uniform  niimber  of  days  each  year  throughout  the  cycle  causes 
the  degree  of  disturbance  represented  by  the  selected  days  to  vary.  Each  method  presents  dif- 
.  ficulties  in  evaluating  the  time  variation  of  the  recurrence  tendency. 


Evaluation  of  the  recurrence  tendency  may  be  accomplished  by  calculating  coefficients  of 
linear  correlation  between  the  magnetic  character  figure  for  each  day  and  the  character  figure 
for  the  nth  succeeding  day.  Thus 
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where  rQ,  n>  the  auto -correlation  coefficient  between  day  0  and  n,  is  derived  from  the  character 
figures  x  expressed  in  units  of  their  standard  deviation  cO*.  The  recurrence  measured  in  this  way 
will  apply  for  all  character  figiures  rather  than  for  some  designated  degree  of  disturbance.  Year 
to  year  comparison  of  the  recurrence  tendency  is  effected  by  calculation  of  the  coefficients 
separately  for  each  year.  Auto -correlation  coefficients  showing  this  recurrence  were  computed 
for  the  International  Magnetic  Character  Figures,  C,  which  are  listed  by  CHAPMAN  and  BARTELS 
[1940]  for  1890  to  1937  and  tabulated  periodically  in  the  Journal  of  Terrestrial  Magnetism  and 
Atxnospheric  Electricity  for  1938  to  1944.  The  index  C  rates  every  day  on  a  scale  of  increasing 
magnetic  disturbance  from  0.0  to  2.0,  a  range  convenient  for  correlation  analyses.  The  coefficients 
were  computed  for  each  year  separately;  thus  each  group  has  a  sufficiently  large  population  that, 
although  the  distribution  of  C -figures  is  not  normal  and  varies  slightly  from  year  to  year,  regroup¬ 
ing  seemed  unnecesary.  The  computations  were  carried  out  by  the  Mendenhall-Warren-Hollerith 
correlation  method,  [WARREN  and  MENDENHALL,  1929],  an  exact  method  employing  punch  card 
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machines.  The  correlation  coefficients  for  1-,  2-,  3-, . . .  36-day  intervals  (with  the  exception  of  18- 
and  22-)  and  for  52-,  53-, . . .  56-  and  79-,  80-, . . .  83-day  intervals  for  each  year  1890  to  1944  are 
listed  in  Table  1.  Negative  coefficients  are  over  scored  and  all  coefficients  for  convenience  are 
multiplied  by  100.  Since  the  correlation  coefficients  were  computed  for  several  Intervals  simul¬ 
taneously,  the  ‘V^ar^'  for  which  recurrence  was  measured  does  not  always  correspond  to  the  calendar 
year,  the  difference  being  as  much  as  35  days  in  addition  to  the  length  of  the  interval.  The  initial 
days  of  the  27-,  54-,  and  8-day  intervals,  however,  comprise  a  calendary  year.  The  errors  thus 
introduced  are  unimportant  since  the  coefficients  will  not  be  treated  separately. 

The  characteristics  of  the  recurrence  tendency  of  magnetic  activity  are  well  demonstrated  by 
these  auto -correlation  coefficients.  Figure  1  shows  the  variation  in  the  strength  of  the  tendency 
during  the  55  years,  represented  by  the.  annual  coefficients  for  the  27-day  interval;  these  are  the 
largest  in  the  first  recurrence  pulse  in  more  than  80  per  cent  of  the  individual  years.  The  cyclic 
variation  of  the  recurrence  tendency  is  more  apparent  if  this  curve  is  smoothed  according  to 
fq  =  (r.i  +  2ro  +  r+i)/4  (where  the  subscripts  refer  to  years)  as  in  Figure  2A.  We  recall  that 

the  division  of  the  data  into  approximately  cal¬ 
endar  years  was  arbitrary.  The  recurrence 
cycle  is  approximately  180  degrees  out  of  phase 
with  the  curve  of  sunspot  numbers  (see  Fig.  2B) 
and  also  with  the  curve  of  the  M^-measure  of  mag¬ 
netic  activity  (Fig.  2C),  but  evidently  closely  re¬ 
lated  to  both.  The  recurrence  cycle  and  the 
Fig.  l--27-day  recurrence  correlation  coeffic-  annual  percentage  of  low  latitude  sunspots  (^10°) 
ients,  1890-1944  compiled  from  the  Greenwich  publications,  are 

remarkably  similar.  The  correlation  between 
the  unsmoothed  27-day  recurrence  and  the  others 
is  r  =  -0.56,  -0.4Q,  and  -(-0.52,  respectively. 

To  study  variation  in  strength  of  recurrence 
during  a  single  solar  cycle,  we  may  consider  the 
average  cycle  to  be  made  up  of  three  years  be¬ 
tween  minimum  and  maximum  year  and  six  be¬ 
tween  maximum  and  minimum,  typical  for  the 
last  eight  sunspot  cycles.  Averaging  the  27-day 
recurrence  coefficients  for  the  maximum  and 
minimum  years  and  distributing  the  coefficients 
of  the  remainder  of  the  years  proportionately 
according  to  this  scheme,  we  obtain  Figure  3, 
which  shows  that  rec^l^rence  is  most  pronounced 
about  a  year  before  minimum.  It  is  generally 
strong  from  about  two  years  after  maximum 
through  the  minimum  year,  arid  generally  weak 
during  the  remainder  of  the  cycle.  For  compari¬ 
son,  the  annual  sunspot  numbers  are  also  shown 
in  Figure  3  for  the  same  period,  combined  in  the 
same  way. 

The  variation  from  year  to  year  of  the 
smoothed  coefficients  for  54-  and  81 -day  intervals 
follows  closely  that  for  the  27-day  interval  (see 
Fig.  4).  The  relative  magnitudes  of  the  0-,  27-,  54-,  and  81-day  recurrence  pulses  for  the  stages 
strong  and  weak  recurrence  are  shown  in  Figure  5. 

Existence  of  any  significant  recurrence  intervals  other  than  about  27-days  or  multiples  thereof 
should  be  detectable  in  the  results  given  in  Table  1.  The  annual  coefficients  for  all  intervals  are 
averaged  in  Figure  6  for  four  periods  of  sunspot  activity  .  .  .  minimum  (calendar  year  of  epoch 
+  one  year),  maximum,  and  the  intervening  ascending  and  descending  years  of  the  sunspot  cycle. 

The  very  strong  maximum  in  the  coefficients  at  about  27  days  is  clearly  shown,  but  no  other  pro¬ 
nounced  maximum  is  revealed.  In  the  descending  and  minimum  periods  of  solar  activity  the  co¬ 
efficients  for  13-,  14-,  and  15-day  intervals  are  slightly  larger  than  on  adjacent  days,  about  zero 
instead  of  slightly  negative.  Reference  to  Table  1  shows  that  in  a  few  individual  years,  for  in¬ 
stance  1895,  1922,  1930,  and  1942,  the  coefficients  were  higher  than  average  at  about  these  inter¬ 
vals.  Possibly  this  illustrates  a  tendency  for  active  solar  regions  to  be  situated  180  degrees  ^ 
apart  in  longitude  at  these  stages  of  the  solar  cycle.  A  few  coefficients  of  distinctly  above  average 


Fig.  2--Comparison  of  data- relative  to  27-day 
recurrence  tendency,  1890-1944 
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value  for  intervals  other  than  those  of  recognized 
recurrence  appear  in  Table  1,  such  as  for  the 
six-day  interval  in  1900,  the  20-day  interval  in 
1911,  and  the  17-day  interval  in  1931.  The  evi¬ 
dence  for  consistent  recurrence  intervals  other  than  those  of  about  27  days  or  multiples  thereof 
is  very  weak;  certainly  no  other  interval  is  a  rival  in  strength  and  prominence. 


Fig,  5 — First,  second,  and  third  reciirrence 
pulses  in  years  of- strong  and  weak  recur¬ 
rence  correlation 


Consideration  of  the  one-,  two-,  and  three-day  recurrence  coefficients  (see  Table  1  and  Fig.  6), 
shows  that  there  is  a  very  high  correlation  between  the  magnetic  activity  of  one  day  and  the  first 
following,  a  much  lower  correlation  with  the  second  and  very  little  with  the  third.  Thus,  the  time 
unit  of  geomagnetic  activity  would  seem  to  be  on  the  average  two  to  three  days.  Because  Intervals 
are  considered,  the  initial  part  of  the  curve  of  Figure  6  is  not  doubled  back  on  Itself  In  order  to 
find  the  average  length  of  a  disturbance.  The  two-  to  three-day  geomagnetic  time  unit  suggests, 
for  instance,  that  forecasts  of  magnetic  activity  should  not  be  made  for  daily  Intervals;  preferably 
^  three-day  average  character  figure  estimates  should  be  attempted. 

Application  to  prediction— The  correlation  coefficients  discussed  in  the  preceding  section  may 
be  applied  to  the  prediction  of  magnetic  character  figures,  in  particular  to  the  prediction  of  mag¬ 
netic  storms  and  the  closely  allied  ionospheric  disturbances.  The  problem  is  to  predict  a  char¬ 
acter  figure,  x'q  =  ax^  +  bx2  +  . . .,  for  the  value  xq,  from  the  character  figures  of  several  pre¬ 
vious  days,  X2  .  .  such  that  the  squares  of  the  errors  of  a  large  niunber  of  predictions, 

I>(xq  -  x'q)2,  Is  a  minimum.  Substituting  and  differentiating  with  respect  to  the  unknown  constants, 
a,  b,  .  .  we  have  as  the  required  conditions 

£(x()Xi  -  axj  -  bxjxg  -  cx^xg  -...)=  0 
Z(xqX2  -  axjXg  -  bx|  -  cxgXg  -...)=  0  etc. 
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If  the  character  figures,  x,  are  measured  in  terms  of  their  standard  deviations,  then  £x|/N  = 

,.  .  .  =  1.  Likewise  SIxqXj/N  is  the  auto-correlation  coefficient,  ro>l»  ^®tween  Xq  and  Xj^;  similarily 
iSxQXg/N  is  rQ,2,  etc.  Given  the  correlation  coefficients,  we  may  solve  the  simultaneous  equations 
for  af  b,  .  .  .,  the  constants  for  the  least  squares  prediction  formula  xj  =  axj  +  bx2  +  .  .  . 

Theoretically  any  number  of  terms  may  be  included  in  a  prediction  formula.  The  accuracy  of 
the  prediction  is  not  necessarily  reduced  by  the  inclusion  of  terms  with  small  coefficients,  or,  in 
other  words,  terms  based  on  weak  correlation.  The  labor  of  solving  a  large  number  of  simul¬ 
taneous  equations,  however,  is  considerable,  and  so  we  will  confine  ourselves  to  formulas  of  five 
terms  or  less,  including  the  more  significant  recurrence  intervals  of  1,  26,  27,  28,  and  54  days. 

The  prediction  method  may  also  be  extended  to  include  nonlinear  terms,  for  instance  ones 
involving  the  consistency  of  the  27-day  recurrence  sequence  of  which  the  predicted  day  is  a  member. 
Thus,  if  the  character  figures  on  both  the  27  and  54  preceding  days  were  small,  one  might  expect 
a  greater  probability  that  the  character  of  the  “zero**  day  will  also  be  small.  Accordingly  dis¬ 
turbance  indices  Cbc  =>tc.27>f  C.54)  and  =V(C_27  x  c5_8i)  were  derived  for  each  day,  and 
rQ  j5and  the  appropriate  cross -correlation  coefficients  were  computed  so  that  terms  involving 
Xf^and  X0  could  be  included  in  prediction  formulas  with  X20,  X27,  and  so  forth. 

Because  of  the  demonstrated  variation  in  the  strength  of  recurrence  with  the  solar  cycle, 
separate  prediction  formulas  are  given  for  the  stages  of  strong  and  weak  recurrence.  The  average 
auto -correlation  coefficients  for  the  significant  recurrence-intervals  and  rQ^a  and  rQ  /3are  given 
for  each  subdivision  of  the  cycle  in  Table  2,  With  these  and  the  required  cross-corr elation  co¬ 
efficients,  prediction  formulas  applicable  to  years  of  strong  (1)  and  weak  (2)  recurrence  are: 


X  Q  -  0.23x  +  0,44x2^  +  0.11x2g  +  0,12x2^  +  (Al) 

x'q  =  0.33x  +  0.47x1  +  0.06x26  +  0.07x27  +  0.03x23  +  0.04xg^ . (A2) 

x'q  =  0.46x  +  0.13x26  +  ^•^'^^27'^  ®*^^X23  +  (Bl) 

x^Q  =  0,71x  +  0.06x26  +  0.09x27  +  0,10x2g  +  0»04Xg^  . . 

x'q  =  0.49x  +  0.14x26  +  0.07x27  +  0,13X23  (Cl) 


Table  2 --Average  auto -correlation  coefficients  for  significant  intervals 
in  years  of  strong  and  weak  recurrence  condition 


Type  of  year 

Recurrence  interval  in  days 

1 

1 

27 

28 

54 

81 

a 

0 

Strong  recurrence 

0.51 

0.27 

0.34 

0.27 

0.24 

0.13 

0.35 

0.28 

Weak  recurrence 

0.49 

0.13 

0.18 

0.16 

0.08 

0.04 

0.18 

0.15 

The  subscripts  refer,  of  course,  to  preceding  days,  and  a  is  as  defined  above.  The  units  of  x  are 
those  of  the  magnetic  character  figures,  since  conversion  from  units  of  standard  deviation  is  ef¬ 
fected  by  the  term  in  x,  the  mean  value  of  the  character  figures,  whose  coefficient  is  1  -  (a  +  b  +  . , 

In  formula  (Cl)  the  effect  of  the  strong  27-day  recurrence  correlation  is  shared  by  X27  and 
x^.  Because  the  non-linear  term  receives  the  greater  weight,  we  may  infer  that  our  supposition 
is  valid  that  the  consistency  of  the  27-day  recurrence  sequence  has  importance  in  prediction.  In  a 
formula  like  (Cl),  but  also  including  a  term  in  x^,  the  coefficient  in  that  term  proved  to  be  slightly 
negative.  The  coefficient  of  x,  however,  was  the  same  as  in  formula  (Cl),  so  we  may  conclude 
that  XA,  though  closely  correlated  with  xq,  was  entirely  dependent  on  X27  and  x^^  and  therefore 
contributed  nothing  of  its  own  to  the  prediction. 

The  usefulness  as  well  as  the  limitations  of  this  method  of  predicting  magnetic  character  fig¬ 
ures  is  illustrated  by  trial  predictions  in  typical  years  of  strong  and  weak  recurrence.  Correlation 
of  observed  character  figures  with  the  predictions  made  by  formula  (Al)  for  1944  gave  r  =  0.56 
(see  Fig.  7),  by  formula  (Bl),  r  =  0,37,  while  for  the  same  period  the  simple  27-day  correlation 
coefficient,  rQ  27  =  Predictions  by  formida  (B2)  for  1945  compared  with  observed  values 

gave  r  =  0.17,  *when  rQ  27  =  0,14.  The  prediction  formulas  work  out  better  than  simple  27-day  pre¬ 
dictions  as  expected,  tfie  improvement  being  obtained  mainly  by  including  in  the  prediction  data 
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from  the  day  first  preceding  the  predicted  one.* 
Regression  towards  the  mean  limits  the  range 
of  the  predicted  character  figures  such  that  the 
predictions  for  years  of  weak  recurrence  lack 
detail.  Independent  correlations  of  some  mag¬ 
nitude  used  in  the  formulas  will  increase  the 
range  and  accuracy  of  the  predictions.  It  might 
be  expected,  for  instance,  that  other  relation¬ 
ships  such  as  solar-geomagnetic  correlations 
could  be  used  for  more  reliable  prediction  for¬ 
mulas. 


Fig.  7- -Correlation  diagram  of  magnetic  char¬ 
acter  figures,  1944,  observed  and  predicted 
by  Formula  Al;  r  =  0.56 
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